Low-lying spectroscopy of a few even-even silicon isotopes investigated by means of 
the multiparticle-multihole Gogny energy density functional 
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A multiconfiguration microscopic method has been apphed with the Gogny effective interaction to 
the calculation of low- lying positive-parity states in even-even ^^~^^Si isotopes. The aim of the study 
is to compare the results of this approach with those of a standard method of GCM type and to get 
insight into the predictive power of multiconfiguration methods employed with effective nucleon- 
nucleon force taylored to mean- field calculations. It is found that the multiconfiguration approach 
leads to an excellent description of the low- lying spectroscopy of ^^Si, ^^Si and ^^Si, but gives a 
systematic energy shift in "^^Si. A careful analysis of this phenomenon shows that this discrepancy 
originates from too large proton-neutron matrix elements supplied by the Gogny interaction at the 
level of the approximate resolution of the multiparticle-multihole configuration mixing method done 
in the present study. These proton-neutron matrix elements enter in the definition of both single- 
particle orbital energies and coupling matrix elements. Finally, a statistical analysis of highly excited 
configurations in ^^Si is performed, revealing exponential convergence in agreement with previous 
work in the context of the shell model approach. This latter result provides strong arguments 
towards an implicit treatment of highly excited configurations. 
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I. INTRODUCTION 

Mean field approaches provide a reliable foundation for 
an approximate solution of the nuclear many-body prob- 
lem. Nowadays, a lot of efforts is applied to move beyond 
the mean field approximation and account for missing 
correlations. Special attention is paid to the restora- 
tion of symmetries broken in mean field approaches [l|-[9| , 
for example using projection techniques. An alternative 
way is to develop a theory in which the trial wave func- 
tions preserve certain symmetries. In particular, this 
can be achieved by multiconfiguration methods widely 
used in various applications, including atomic, molecu- 
lar, and condensed matter physics. When the interaction 
is known, this kind of approach provides a ver y a ccu- 
rate description of a system. In a previous work we 
proposed in the nuclear physics context a variational ap- 
proach based on multiparticle-multihole {mp — mh) pro- 
ton and neutron configuration mixings that uses the two- 
body finite-range density-dependent Gogny interaction in 
a fully microscopic way. As a first step, the method was 
applied in a particular case: pairing correlations (proton- 
proton and neutron-neutron) were investigated for the 
ground states of several even-even tin isotopes. A pio- 
neering work using the Skyrme SIII interaction for the 
mean field part and a schematic contact interaction for 
the residual part has been directed to the description of 
i^^- isomers in the ^^^Hf mass region |12l] . 

The main objective of the present study is twofold: 
first, to test the ability of mp — mh multiconfiguration 
approach to describe low-lying nuclear states and second, 
to discuss the statistical features of highly excited con- 
figurations in nuclei. With this aim in view, a few silicon 
isotopes have been chosen and their excited states cal- 



culated with the mp — mh approach. Concerning highly 
excited configurations, the leading idea is to investigate 
whether exponential convergence behavior, revealed in 
the standard shell model (SM) approach p!3l-[T7|, also 
emerges from variational {mp — mh) configuration mixing 
methods in which single particle excitations are not 
restricted to a single major shell. 

As practical calculations inevitably require some trun- 
cation of the orbital space and order of excitation, the 
mp—mh method proposes a promising scheme that allows 
one to predict the energies of low-lying states in a very 
accurate way. The possible use of statistical properties 
of highly excited states, which display generic signatures 
of quantum chaos close to random matrix theory, dras- 
tically reduces the sizes of the explicit diagonalizations. 
In the literature, one finds other approaches proposing 
different algorithms aSjfor example, the density matrix 
renormalization group |19l-[23| or Monte Carlo techniques 
[2^ , selecting the most relevant configurations for the de- 
scription of many-body states. 

The paper is organized as follows. The main charac- 
teristics of the mp — mh multiconfiguration approach ap- 
plied in this study are presented in Section [TTl In Section 
Hm the results obtained with the mp — mh approach for 
low- lying states in ^^~^^Si are presented and compared 
with those derived from a five- dimensional approximate 
Generator Coordinate Method (GCM) approach (subsec- 
tions [IIlAl and [IirB]). Differences between theoretical and 
experimental results are discussed. In particular, the sys- 
tematic energy shift found in ^^Si is analyzed in terms of 
proton-neutron matrix elements. Subsection IIII CI high- 
lights the crucial role played by the residual interaction 
between protons and neutrons. Section HVl is devoted to 
the analysis of the statistical properties of highly excited 
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configurations, taking ^^Si as an example. Conclusions 
and perspectives are given in Section IVl 

II. MULTIPARTICLE-MULTIHOLE 
CONFIGURATION MIXING APPROACH 

In this part, we discuss a few characteristics of the 
mp — mh configuration mixing approach. A general 
derivation of the corresponding formalism in nuclear 
physics, with two-body density-dependent interactions, 
has been introduced in Ref. [nj. It is worth to recall 
here the basics of the method not only to fix notations 
but also to provide an alternative formulation of equa- 
tions in terms of a "core + valence space" description. 

In the mp — mh configuration mixing method, the ef- 
fective Hamiltonian is defined as a functional 

H{p) = k + V{p) (1) 

of the single-particle density matrix p. In Eq. ([1]), the 
Hamiltonian contains a kinetic term K (which includes 
the one-body center-of-mass corrections) and a two-body 
density-dependent potential term V{p) (which includes 
the Coulomb potential for protons as well as the two- 
body center-of-mass corrections). In our study, the DIS 
Gogny interaction [l8| is adopted. 

The trial wave functions |^) that describe nuclear sta- 
tionary states are expressed as linear combinations 

1^) = ^c.^c.Mc.Ac..) (2) 

of direct products 

l^a^^aj = l^a^) ^ |(/>aj (3) 

of proton and neutron Slater determinants, |(/>q;^) and 
|0Q,^) respectively, containing a priori any multiple p-h 
excitations that respect conserved quantum numbers. 

Eq. ([2]) contains two sets of unknown parameters, the 
mixing coefficients A^^^^ and the single-particle orbitals. 
Both are determined by applying a variational principle 
to a functional T{p) related to the total energy of the 
system, 

T{p) = imipm - A{*i*> - E (4) 

where A and are Lagrange multipliers and Qi possi- 
ble additional constraints that we will leave out in the 
following. One assumes that the one-body density p en- 
tering the effective Hamiltonian H{p) is the correlated 
one: p = 

The minimization of J^{p) with respect to the A^^ar, 
leads to a non-linear secular equation that is equivalent 
to a diagonalization problem in the multiconfigurational 
space of a Hamiltonian matrix 

^ ^ "^a-jj^ai, ,a'^a'^AQ,'^Q,'^ = XAq,^q,^ . (5) 



In Eq. (j5j), the matrix H contains contributions of the 
Hamiltonian H{p) and of rearrangement terms that come 
from the density dependence of the interaction, 

mnr 

(6) 

where the summation over r specifies the proton and 
neutron contributions to the generalized rearrangement 
terms with coefficients 71^^. It is the presence of the 
rearrangement terms that transforms Eq. ([5]) into a 
non-linear eigenvalue problem. As can be seen from Eq. 
()A5p . the rearrangement terms contain contributions as- 
sociated with the one-body density p and the two-body 
correlation matrix a defined by 

cTijM = (^|a+a^a/a^ |^) - pjipik ^ pjkpii- (7) 

In Appendix [XI we express Eqs. ([5])- ([6]) in a "core + 
valence space" scheme, a form which is explicitly used in 
the present study. 

The minimization of J^{p) with respect to the single- 
particle orbitals leads to inhomogeneous Hartree-Fock 
(HE) equations which depend on the amount of corre- 
lations contained in |^) 

[h{p,a),p] = G{a), (8) 

where 

Gkiicr) = ^ ^{im\V{p)\kn)aii^mn 

1 ^ (9) 

imn 

In Eq.(|8]), h{p^a) is the one-body mean-field Hamilto- 
nian built with the one-body density p and the two-body 
correlation matrix a: 

h,j{p,a) = {i\K\j) ^T,j{p) ^ dT,j{p) ^ dT,j{a). (10) 

Explicit expressions for the fields rij(p), dTij{p) and 
dTij{a) are given in Ref. The basis that diagonal- 

izes h{p^ a) provides proton and neutron single-particle 
orbitals with energies e[. The diagonal part of Eq. (j6j), 
that is the one obtained by taking a = a' ^ can be easily 
written in terms of the ej . In the limit where the mixing 
reduces to the sole HF configuration, one recovers the 
standard HF expression for e[, 

, (11) 

t't" h' = l h" — l 

In Eq. (pTj) . summations over /i, h' and h" run over the 
hole states. 
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FIG. 1: (Color online) Triaxial HFB potential energy surfaces for ^^Si (panels (a), (b), (c) and (d), respectively) obtained 
with the DIS Gogny interaction. 



As discussed in Ref. fH] , a fully self-consistent solution 
of the multiconfiguration approach consists of solving si- 
multaneously both equations (jSj) and ([8]). It is important 
to notice that, in our approach, the mean field and be- 
yond mean-field descriptions are obtained in a consistent 
way since they both follow from variations of the same 
functional. A standard way to solve both equations is 
to use an iterative procedure. Starting from a HF so- 
lution (characterized by p = phf and a = 0) that pro- 
vides an initial set of single particle orbitals, one builds 
the multiconfiguration basis and solves the configuration 
mixing Eq. (|5]). Then, from the mixing coefficients ob- 
tained from Eq. (|5j), one calculates p and a and solves 
Eq. dH]). With the new set of single particle orbitals, one 
redoes the procedure until convergence of p and a. In 
this way, the single-particle orbitals contain effects com- 
ing not only from the mean-field built with the one-body 
density matrix p but also from the correlations contained 
in the matrix a. 

In general, the actual solution of Eq. (j8j) is a very diffi- 
cult task. In a previous study dedicated to the particular 
case of pairing-type correlations in the ground states of 
even-even tin isotopes one has solved Eq. (|8|) in an 
approximate way, by neglecting a. In the present study of 
even-even silicon isotopes, as a first step, we have solved 
only the configuration interaction (CI) part (Eq. (j5j)) of 
the multiconfiguration approach. 



As already pointed out in Introduction, multiconfigu- 
ration methods are able to preserve certain symmetries. 
Concerning particle numbers, the correlated wave func- 
tion ([2]) is a superposition of Slater determinants that 
conserves exactly the numbers of protons and neutrons. 
As to angular momentum, the multiconfiguration equa- 
tions have been solved in the present study only for spher- 
ical nuclear configurations, with valence spaces compris- 
ing complete spherical subshells. As a consequence, nu- 
clear states have a good total angular momentum J. Ac- 
tually, since our formalism has been developed in axial 
symmetry in order to be able to introduce quadrupole 
deformation explicitly, the only conserved quantum num- 
bers we have are the projection K of the total angular 
momentum and the parity. The conserved quantum num- 
ber J is then obtained in the usual manner by perform- 
ing calculations for successive values K = 0, 1, 2, . . . and 
identifying degenerate multiplets. Let us note that the 
finite axial harmonic oscillator bases used for expand- 
ing single-particle states (see below) are defined with a 
spherical truncation. 

III. LOW-LYING STATES IN 26-32 gj 

This section is devoted to the description of the low- 
lying spectroscopy of ^^~^^Si using the mp — mh con- 
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figuration mixing approach. In the first part of this 
analysis, we investigate the mean- field properties of the 
ground states provided by the Hartree-Fock-Bogoliubov 
(HFB) approach using the same DIS Gogny interaction 
and compare them with the results of a five-dimensional 
approximate GCM approach. 

Technically, the single-particle states introduced in the 
different approaches (HF, HFB, mp — mh) are expanded 
onto the harmonic oscillator (HO) bases. In the present 
work, 7Vo=ll major spherical shells have been taken. 
This basis size has been found sufficient to ensure the 
convergence of all the results obtained in this work in 
the three approaches mentioned above. For instance, the 
convergence of low-lying state energies in the mp — mh 
configuration mixing approach has been achieved within 
an accuracy better than 0.1 keV. 



A. Ground state deformation properties 

In order to investigate the properties of the ground 
states of ^^~^^Si, we start with triaxial HFB calculations 
constrained according to the dimensionless deformation 
parameters f3 and 7, 
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and 7 = tan 



(12) 



r), q2 = (Q2) 



In Eq. (Ha), go = (Qo) = 
(x^ — y'^) and ro=1.2 fm. 

Fig. [1] displays the triaxial HFB potential energy sur- 
faces (PES) of ^^~^^Si in the P and 7 degrees of freedom. 
Isolines associated with total energy are indicated with a 
numbering corresponding to the height of the potential 
(in MeV) relative to the minimum of the HFB potential 
(dark blue online) for each nucleus. The /3hfb and 7hfb 
values for the HFB solution of lowest energy are indicated 
in the second and the third columns of Table HI 



Nucleus /3hfb 


7HFB 


(/3)5DCH 


(7)5DCH 


2^Si 


0.32 


60.0° 


0.41 


28° 


''Si 


0.37 


60.0° 


0.40 


27° 


3°Si 


0.00 


0.0° 


0.39 


29° 


^^Si 


0.01 


34.0° 


0.37 


28° 



TABLE I: /3 and 7 deformation properties of the ^^~^^Si 
ground states from the HFB and 5DCH approaches and the 
DIS Gogny interaction. 

From Table IH one sees that the HFB ground states 
of ^^Si and ^^Si are similarly characterized by a large 
value of /3hfb — 0.35 and 7hfb = 60°. The HFB minima 
are found well-deformed on the oblate side. ^^Si and ^^Si 
exhibit different characteristics. Pufb is equal to zero for 
^^Si and very close to zero for ^^Si. The values of 7hfb 
are very different, 7hfb = 0° for ^^Si and 7hfb = 34° 



for ^^Si. They can be considered as spherical and nearly 
spherical nuclei, respectively. 

Investigating the PESs of Fig. (TJ one can deduce that, 
even though /3hfb and 7hfb may be quite different in 
the four nuclei, the common feature of the four PESs is 
their softness in both the /3 and 7 degrees of freedom. In 
relative, the ^^Si PES appears to be more 7— soft than 
the one of ^^Si. The PESs of ^^Si displays similar fea- 
tures as the one of ^^Si. This proposes an important role 
of triaxiality. This softness is quantitatively evidenced 
by the mean values (/3)5dch and (7)5dch displayed in 
Table HI These values have been obtained from a five- 
dimensional collective Hamiltonian (5DCH) describing 
both /3 — 7 and rotation modes with the use of the ap- 
proach developed in Ref. [ol, [lo| . Let us recall that such 
an approach is based on the completely microscopic gen- 
erator coordinate method (GCM) [l| and allows one to 
find collective excitations of pure rotational-vibrational 
character from the only data of the nucleon-nucleon ef- 
fective force. The 5DCH calculations predict strong dy- 
namical /3-deformations for the ground states of all four 
silicon isotopes together with significant triaxiality. One 
observes that the 5DCH collective dynamics introduces 
considerable changes with respect to HFB in the /3 — 7 
ground state deformations. 

To detail the information provided by (/3)5dch and 
(7)5DCH mean values. Fig. [2] shows the 5DCH collective 
wave functions of ^^~^^Si ground states (panels (a), (b), 
(c) and (d), respectively). Even though the ^^Si collective 
wave function is partly suppressed for 7 between 45 and 
60 degrees, the spreading of collective wave functions in 
/3 — 7 plane is evident. The spreading in /3 goes up to 
-0.8 for 26Si and ^^Si and -0.6 for ^^Si and ^^Si. 





7rd5/2 


vrsi/2 


7rd3/2 


Z^Pl/2 


Z^d5/2 


Z^Si/2 


Z^d3/2 


Z^f7/2 


2^Si 


-7.07 


-2.61 


0.88 


-27.73 


-15.96 


-10.73 


-7.33 


-1.37 


''Si 


-10.05 


-5.16 


-1.69 


-28.39 


-15.95 


-10.91 


-7.34 


-1.18 


3°Si 


-12.60 


-7.32 


-4.11 


-26.75 


-16.21 


-11.56 


-7.29 


-1.41 


^^Si 


-15.09 


-10.08 


-7.51 


-26.44 


-16.28 


-11.40 


-7.80 


-2.01 



TABLE II: Energies (in MeV) of spherical proton and neutron 
HFB single-particle states in ^^~^^Si. 

In Table |TT1 we report the energies ej of the spherical 
orbitals for ^^~^^Si, for protons {7rld^/2^ 7r2si/2, '^^^(^3/2) 
and neutrons (i^lpi/2, ^l<^5/2, ^2^1/2, ^^1/7/2)- 
These energies slowly evolve from one isotope to the 
other. We notice a slight increase of the neutron gap be- 
tween z^2si/2 and orbitals in ^^Si, which reaches 
a maximum value of 4.56 MeV. A previous study of 
TV = 16 isotones with the DIS Gogny interaction have 
suggested that all Z = 10 — 18 isotones show strong de- 
formations, limiting the understanding of A/" = 16 as a 
magic number to the sole oxygen neutron drip line ^5| . 

The evolution of proton and neutron single-particle or- 
bitals obtained within the HFB approximation with axial 
deformation /3 (7 = 0°) is displayed in Fig. [3l To dis- 
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FIG. 2: (Color online) Ground state 5DCH collective wave functions for ^^Si (panels (a), (b), (c) and (d), respectively). 



cuss the general characteristics of single-particle orbitals 
when spherical symmetry is broken, we have arbitrarily 
selected the isotopes ^^Si (for protons) and ^^Si (for neu- 
trons). Open circles (black online), squares (red onhne), 
stars (green online) and triangles (blue online) stand for 
the projections of the angular momentum on the sym- 
metry axis with values jz = 1/2, 3/2, 5/2 and 7/2, re- 
spectively. Solid (dashed) lines correspond to positive 
(negative) parity orbitals. The chemical potential is in- 
dicated with filled black circles and is denoted by A. 

For protons (upper panel), starting from the oblate 
side and up to the spherical shape, the Fermi level corre- 
sponds to the jz=l/2 deformed orbital originating from 
the shell. From sphericity up to a prolate defor- 

mation P :^0.5, the Fermi level is located on the jz=^/2 
deformed orbital and then migrates to the jz=l/'2 orbital 
coming from the 2si/2 shell; finally, for a very large value 
of /3, it follows the jz=l/2 deformed orbital coming from 
the 1/7/2 shell. The ^^=1/2 and jz=3/2 deformed or- 
bitals from the shell for small oblate deformations 
fall down instead of going up. Investigating this plot in 
more detail, one sees that similar trends are encountered 
for deformed orbitals originated from the 2si/2 and lds/2 
shells. One is led to the conclusion that these orbitals are 
strongly mixed through deformation as confirmed by the 
presence of avoided level crossings in single-particle spec- 



tra. The natural continuity of the oblate ^^=1/2 from 
the 1(^5/2 shell is the prolate ^^=1/2 from the 1(^3/2 shell. 
Other examples are the oblate ^^=1/2 from the 2^1/2 
shell and the prolate ^^=1/2 from the shell, the 

oblate jz=l/2 level from the lds/2 shell and the prolate 
jz=l/'2 from the 2si/2 shell, or the oblate/prolate jz=3/2 
from the ld^/2 shell and the prolate/oblate jz=3/2 from 
the ld^/2 shell. From this kind of analysis, one ex- 
pects that the proton excitations contributing to low- 
lying states will mainly arise from the sd-shell, while the 
influence of the ^^=1/2 level from the I/7/2 shell will ap- 
pear only at quite large P deformation. It is important 
to note that these mixings and repulsions of deformed 
orbitals are not necessarily synonymous with inversions 
of orbitals. 



The behavior of the neutron orbitals is similar. As in 
the present paper we study isotopes where the 2si/2 and 
1(^3/2 neutron subshells are filled or partially filled, it is 
important to look at higher shells. Here we only men- 
tion that the same game as for protons is played between 
deformed neutron orbitals originating from I/7/2, 2^3/2 
and 2pi/2 with low jz values, see Fig. [3l lower panel, and 
that upper shells of negative parity play a role at small 
deformations only for the heavier isotopes. 
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FIG. 3: (Color online) Evolution of proton and neutron single- 
particle orbitals obtained within the HFB approximation with 
the axial deformation /3 (7 = 0°) in ^^Si (upper panel (a)) and 
^2 Si (lower panel (b)). 



B. Low-lying spectroscopy 



In this subsection we discuss the low-lying states in the 
even-even ^^~^^Si isotopes as predicted by the mp — mh 
configuration mixing approach and compare them with 
the ones obtained with the 5DCH method. Only posi- 
tive parity states are investigated in this work. Exper- 
imentally, these nuclei are challenging as they exhibit a 
large variety of states at low energy, and strong changes 
arise from one isotope to another (for example, in the 
state). For this reason, from the theoretical view- 



point, they can be considered as benchmarks for both the 
many-body method employed and the properties of the 
effective nucleon-nucleon interaction used. The compari- 
son with the 5DCH approach, a method that has proved 
its pertinence through a global survey [9], is interesting 
in the sense that, as the same DIS Gogny interaction 
is used, it enables one to specify the role of rotational 
and quadrupole correlations in the spectroscopy of the 
positive parity states in silicon isotopes. 

All the mp—mh results we present below correspond to 
mixing within the sd-shell. We have checked the influence 
on the low-lying states of our interest when including the 
Ip and If subshells. In particular, the adding of lpi/2 
and lp3/2 had essentially no effect on the spectroscopy 
of ^^Si, ^^Si, ^^Si and ^^Si. For the selected spectroscopy 
of ^^Si, the effect of the I/7/2 shell was very small. For 
^^Si, only small variations were observed and they con- 
cerned only the highest states. The largest differences (~ 
500 keV) due to the I/7/2 and 2^3/2 orbitals have been 
encountered in ^^Si and ^^Si. Concerning the 5DCH ap- 

Nucleus States Experiment mp — mh 5DCH 



2^ 


1.795 


1.502 


2.426 


2J 


2.783 


2.567 


5.124 


ot 


3.332 


3.740 


8.146 




3.756 


3.233 


9.126 


4+ 


3.842 


3.293 


6.119 


3+ 


4.093 


3.779 




2t 


4.138 


3.915 




4+ 


4.183 


4.758 


9.254 


2+ 


4.446 


4.783 




O3+ 


4.806 


4.959 




4+ 


5.229 


4.931 






5.330 


5.741 




2+ 


5.562 


5.577 




3+ 


? 


5.755 




Ot 


5.940 


6.690 




2| 


6.350 


6.823 





2+ 



TABLE III: Excitation energies (in MeV) of positive parity 
low-lying states in 26Si from experiments (second column), 
and calculated with the variational mp — mh configuration 
mixing method using the DIS Gogny force (third column). 
In the fourth column a few energies derived from the 5DCH 
approach are displayed. The symbol * means that the 
spin of the state is not experimentally assigned. The symbol 
? indicates a state which is not seen experimentally. 

proach, only theoretical results with energies lower than 
c^l2 MeV are presented. This choice is not fully arbi- 
trary as it is motivated by the relative maximum height of 
PESs {c^ 28 MeV above the HFB minimum in the present 
case, see Fig. [1]) used to perform the GCM configuration 
mixing of nuclear shapes. 
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9.479 
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9.496 


9.756 





TABLE IV: Same as Table [ml for ^^Si; The symbol * means 
that spin-parity quantum numbers are not assigned experi- 
mentally. 



For the four isotopes, 3 is the lowest negative parity 
state. Its experimental energy is 6.789 MeV in ^^Si, 6.879 
MeV in ^^Si, 5.487 MeV in ^^Si and 5.288 MeV in ^^Si. 
The decrease of its value gives a flavor of the increasing 
importance of negative parity orbitals at low energy in 
^^Si and ^^Si. However, for the description of positive 
parity states, negative parity orbitals should play a role 
at energies higher than for the description of negative 
parity states as one has to introduce at least 2p — 2h 
excitations to produce a positive parity state. 

Excitation energies calculated with the mp — mh and 
5DCH methods for ^^~^^Si are compared to experimen- 
tal values in Tables IlllilVIl (the energies are expressed 
in MeV). In Table IVl a sixth column named mp — mhg 
has been added, showing mp — mh results shifted by 2.5 
MeV. Positive parity has been assumed for the states in 
^^Si whose spins are not assigned experimentally, based 
on the plausible hypothesis that the lowest negative par- 
ity state is the observed 3~ one. The symbol ? means 
that the state is not observed experimentally. For ^^Si, 
the symbol * indicates that both spin and parity have 
not been measured experimentally; again, making the 
assumption of positive parity, the spin is given by our 



model. For ^^Si, because of the energy of the lowest ob- 
served 3~ state, positive parity is assumed for the fourth 
state at 5.220 MeV. Experimentally, the spin is expected 
to be in the range from 1 to 4; our model predicts spin 3. 
Concerning the 4+ state, the experimental assignment is 
either 4+ or 5~. 



Nucleus States Experiment mp — mh 5DCH mp — mhs 



27 


2.235 


4.609 


2.222 


2.109 


2j 


3.498 


5.704 


4.729 


3.208 


1^ 


3.769 


6.338 




3.838 


n+ 


Q TOO 


'7 TOO 

( . ( oz 


r .oil) 


K OQO 


23+ 


4.810 


8.230 


11.832 


5.730 


3^ 


4.831 


6.709 


8.056 


4.209 


3J 


5.231 


7.904 




5.404 


4^ 


5.279 


7.539 


5.691 


5.039 


O3+ 


5.372 


8.950 


8.585 


6.450 


24+ 


5.614 


9.262 




6.762 




5.950 


8.911 


8.714 


6.441 


25+ 


6.538 


10.186 




7.686 


04+ 


6.642 


9.030 




6.530 




6.865 


8.831 




6.331 




6.914 


10.594 




8.094 


5^ 


6.998 


10.926 




8.426 



TABLE V: Same as Table [ml for ^^Si. 

Theoretical results are provided for 20 states in each 
isotope, except for ^^Si where the spin and parity of some 
states with excitation energies larger than ~5.5 MeV 
have not been firmly assigned experimentally. All excited 
configurations from the sd-shell have been introduced in 
the mp — mh wave functions, up to 6p — 6/1 on the proton 
side and 2p — 2/i, Ap — 4/i, Qp — 6h on the neutron side 
depending on the isotopes. Then, the presented results 
contained up to lOp—lOh configurations in ^^Si and ^^Si, 
12p — 12h configurations in ^^Si, and Sp — Sh configura- 
tions in ^^Si. 

In the case of ^^Si, one sees that there is a very good 
agreement between experimental and mp — mh configu- 
ration mixing energies, whatever the spin and the excita- 
tion energy. Concerning the 5DCH approach, the energy 
of the 2^ state is found too high by ~ 700 keV and the 
energies of the other excited states are also strongly over- 
estimated (by several MeV). In the case of ^^Si, despite 
the inversion between the 4f and O2" levels, quite good 
agreement with experiment is obtained by the mp — mh 
approach. Again, the 5DCH approach tends to overes- 
timate excitation energies. However, we note that both 
theoretical approaches describe the 2^ state within the 
same accuracy as in ^^Si and that the energy of the 22 
state is particularly well reproduced. 

In ^^Si, as seen from Table [V| the energy of the 2^ 
is overestimated by ~2.5 MeV with the mp — mh con- 
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figuration mixing. What is even more surprising is tiiat 
all states appear to be shifted upwards. Actually, reduc- 
tion of all the excitation energies by ~2.5 MeV gives a 
much better agreement with experiment. This can be 
seen from the last column "mp — m/ig" of Table IVl that 
gives the values provided by the mp — mh configuration 
mixing approach minus 2.5 MeV. Then, the discrepancy 
with experiment is reduced to ~0.69 MeV (on average), 
and the theoretical level sequence becomes very similar to 
the experimental one. We have checked that this global 
shift cannot be removed by adding the I/7/2 shell to the 
valence space. Its origin will be discussed later in this 
subsection. 

On the other hand, the energies of several states in this 
nucleus are well reproduced by the 5DCH approach, in 
particular those of the 2^ and 4^ states. Let us note in 
this respect that the 5DCH method does not explicitly 
make use of the matrix elements of the residual inter- 
action between excited configurations, exploiting instead 
the quadrupole deformation properties of the mean field. 
As seen and discussed in relation to Fig. [3l following the 
chemical potential, the deformation is able to catch part 
of the correlation information coming from upper spher- 
ical shells directly or through the mixing of deformed 
orbitals (see discussion in section UlI Al on the evolution 
of single-particle orbitals). In the case of ^^Si, both theo- 
retical methods provide a good description of the selected 
experimental states (except for the O2" with the mp — mh 
configuration mixing method and the 2^ with the 5DCII 
approach) . 

The results obtained with the mp — mh configuration 
mixing indicate that this approach is capable of repro- 
ducing quite well the low-energy spectroscopy of ^^Si and 
^^Si and to a lesser extent, the one of ^^Si. The mean 



Nucleus States Experiment mp — mh 5DCH 



2r 


1.941 


1.959 


2.215 




4.230 


4.871 


5.014 


ot 


4.984 


6.810 


5.318 


3^ 


5.220 


6.004 




2+ 


5.412 


6.758 


9.335 


4+ 


5.502 


6.567 


5.470 



TABLE VI: Same as TableHIIlfor ^^Si. The symbol * indicates 
that spin and/or parity is/are not assigned experimentally. 

deviations between theory (mp — mh) energies and exper- 
imental ones are found to be ~ 369 keV in ^^Si, ~ 653 keV 
in ^^Si and ~ 946 keV in ^^Si. The increase of the devia- 
tion between ^^Si and ^^Si can be attributed to the sys- 
tematic overestimation obtained in the calculation of ^^Si 
highest levels. The same phenomenon holds for ^^Si. As 
pointed out previously, this effect can be partly ascribed 
to the absence of the I/7/2 shell, and, to a smaller extent, 
of the 2p3/2 shell in the present calculations. Nonethe- 



less, the agreement with experiment of mp — mh energies 
can be considered as rather encouraging, considering the 
fact that the DIS Gogny interaction has not been devised 
to describe the kind of general correlations introduced in 
a multiconfiguration approach. In particular, the proton- 
neutron matrix elements between excited configurations 
given by this interaction have not been constrained. 

The following discussion is dedicated to the under- 
standing of the origin of the shift obtained with the 
mp — mh configuration mixing method for excited states 
in ^^Si. In a schematic way, both energy gaps be- 
tween single-particle orbitals and coupling matrix ele- 
ments (ME) between configurations are the key quanti- 
ties that drive the low-lying spectroscopy. One can infer 
that a downward shift can be obtained either by decreas- 
ing gaps (an effect similar to the monopole shifts pointed 
out in the shell model approach [24-26]), and/or by vary- 
ing coupling matrix elements. Proton- neutron matrix el- 
ements are suspected to be mainly responsible for the en- 
ergy shift encountered in ^^Si. In fact, by changing "by 
hand" the values of selected proton-neutron ME implying 
the spin-orbit partners and and using them 

in realistic mp — mh configuration mixing calculations 
(which produces a modification of the energy of the lds/2 
shell), one can derive excited states in ^^Si which are in 
reasonably good agreement with experiment, with devi- 
ations similar to the ones found in ^^Si. It is important 
to note that such changes in matrix elements essentially 
effect only the ^^Si spectrum. In particular no significant 
modification of ^^Si spectrum is observed. In addition, as 
we will see in section HVl where the chaotic behaviour of 
highly excited Slater determinants is studied, too strong 
couplings are found essentially in K = cases, where 
common proton-neutron matrix elements are involved. 

At this stage of our analysis, one has to recall that 
the general formalism of the mp — mh configuration mix- 
ing method exposed in section [III implies that not only 
the secular equation ([5j) has to be solved but also Eq. 
([8|). These two equations in principle provide the "best" 
single-particle representation, i.e. the one that minimizes 
the total energy consistent with correlations. Clearly, the 
solution of Eq. (|8]) may introduce modifications on both 
single-particle energies and coupling ME. However, dis- 
cussing the kind of renormalization produced by Eq. ([8]) 
is far beyond the scope of the present paper, and it will 
be left to a dedicated study. Let us simply mention here 
that, in the context of the present work, introducing a 
tensor term in the effective interaction we would prob- 
ably reproduce the right energy evolution of spin-orbit 
partners In addition, a crude comparison between 
sd-shell ME calculated from the DIS Gogny interaction 
and from the USD interaction used in the shell model 
[3Q| displays large discrepancies essentially in the T = 
channel, where the renormalization effects are expected 
to be the largest. The average difference is equal to ~0.3 
MeV in T = 1 channel and ~1.5 MeV in T = one. 
The large difference found in T = channel is attributed 
essentially to two ME and comes from the lack of ten- 
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sor term by comparing the different contributions to ME 
(central, spin-orbit and tensor). 
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FIG. 4: (Color online) Experimental and calculated {mp — mh 
and 5DCH) E{4:+)/E{2+) ratios. For ^°Si, the shifted for 
mp — mh conflguration mixing energies have been used (see 
text). 



Nucleus 


£(4+) 

E(2+) 
^ 1 ' exp. 


E(4+) 

^ 1 ^ mp — mn 


E(4+) 

^^^1 ) 5DCH 


26Si 


2.140 


2.190 


2.520 


2«Si 


2.595 


2.695 


2.611 


30si 


2.362 


2.389" 


2.561 


32Si 


2.834 


3.350 


2.469 



the experimental trend. For ^°Si, the shifted values for 
the mp — mh configuration mixing approach have been 
used. The significant overestimation in ^^Si, where the 
2f energy is well reproduced, comes from the fact that 
the A'l energy is slightly overestimated. 



C. Role of the proton- neutron residual interaction 

This part discusses the role of the proton-neutron 
residual interaction and the importance of using an ef- 
fective nucleon-nucleon interaction that manifests good 
properties in the T = channel when using mp — mh 
configuration mixing methods. As a benchmark nucleus, 
we have chosen ^^Si. Table IVlnl lists the first seven ex- 
cited states of ^^Si. The experimental values of excitation 
energies and the theoretical ones {mp — mh) are given in 
columns 2 and 3. Columns 4 and 5 display the results ob- 
tained within the mp—mh configuration mixing approach 
when the residual proton-neutron interaction is turned 
off. As the Hamiltonian in use is not exactly isospin- 
invariant (a Coulomb term is included), this symmetry 
breaking leads to a small difference (less than ~ 300 keV) 
for proton and neutron solutions, noted mp — mh^ and 
mp — m/ijy, respectively. 

States Experiment mp — mh mp — mh-K mp — mhu 



2^ 


1.779 


1.993 


5.733 


5.831 




4.618 


5.372 


6.553 


6.712 


ot 


4.980 


4.409 


9.588 


9.651 


St 


6.276 


6.365 


9.732 


? 




6.690 


8.759 


9.873 


9.893 




6.888 


7.769 


? 


? 




7.380 


7.280 


10.283 


10.415 



TABLE VII: Numerical values of experimental and calculated 
{mp - mh and 5DCH) ^(4^)/^(2^) ratios. The superscript 
"s" for mp — mh ratios in ^°Si means that the energies shifted 
by 2.5 MeV have been used, see text. 

In conclusion of this part, one can say that, in the four 
silicon isotopes, most states contain more complex corre- 
lations than the usual collective quadrupole/rotational 
ones. In order to identify the nature of the ground 
state band, we have calculated the energy ratio (4^/2^) 
whose value is a standard indicator of the vibrational, 
rotational or 7-soft nature of nuclei. Results for exper- 
iment, mp — mh configuration mixing and 5DCH meth- 
ods are displayed in Fig. HI Numerical values are re- 
ported in Table IVIII for experiment, mp — mh configu- 
ration mixing, and 5DCH methods. From experiment, 
one observes that ^^Si is close to the vibrational limit 
((4+/2+)vib = 2), while ^^Si and ^^Si are more 7-soft 
nuclei ((4^/2^)^ = 2.5)). In ^^Si, this ratio increases 
towards the rotor limit ((4]^/2]^)rot = 3.3)). The results 
obtained with the mp — mh configuration mixing show 
that data in ^^Si and ^^Si are well reproduced and display 



TABLE VIII: Excitation energies (in MeV) of low- lying states 
in ^^Si calculated with the variational mp — mh configuration 
mixing method with and without residual proton-neutron in- 
teraction (see text). 

The sensitivity of the excitation energies to the proton- 
neutron residual interaction depends on the nature of 
states. For example, this interaction brings the energy 
of the 2t state from ~ 6 MeV to ^ 2 MeV; its impor- 
tance for the structure of correlated wave functions is il- 
lustrated in Table IIXI where the components of the wave 
functions of the and O2" states are listed in two cases. 
We define the quantity Wn that measures the correlation 
content of the wave functions in terms of the order of 
excitation n, namely Op — O/i, Ip — l/i, 2p — 2/i, ... For a 
given eigenfunction l^/?), 

w^ = Y.\4S^ (13) 

where kn represents the Slater determinant components 
with np — nh excitations (n = n^^ + n^y). Case (a) in Ta- 
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State Case Wq Wi W2 Wq W7 

0+ (a) 34.48 0.00 34.77 11.55 12.78 4.25 1.76 0.34 

(b) 93.77 0.00 6.03 0.02 0.17 0.00 0.00 0.00 

0+ (a) 43.70 0.00 12.26 14.12 16.37 8.68 3.82 0.90 

(b^) 0.00 91.51 2.82 5.34 0.20 0.12 0.00 0.00 

(b.) 0.00 91.39 3.25 5.04 0.21 0.11 0.00 0.00 



TABLE IX: Weights Wn (n<8) calculated for 0+ and 0+ 
states in ^^Si. Case (a) corresponds to the full mp — mh calcu- 
lation and case (b) to a calculation without residual proton- 
neutron interaction. The index tt (ly) specifies the proton 
(neutron) solution. 



ble IIXI corresponds to a full calculation and case (b) to a 
calculation without residual proton-neutron interaction. 
The wave functions 0^ and O2" have quite different struc- 
tures. For the ground state, excited 2p—2h configurations 
play a role as important as the initial Op — Oh configu- 
ration. When the proton-neutron residual interaction is 
turned off, most of the correlations disappear. One sees 
that the 2p — 2h configurations built from Ip—lh proton 
excitations combined with Ip—lh neutron excitations are 
essential for the description of the ground state. For the 
excited state, the absence of proton-neutron interaction 
is even worse as it destroys fully the Op — Oh component 
and produces the solutions based on a Ip—lh proton 
or on a Ip — l/i neutron configuration. One observes 
that the proton-neutron interaction brings a lot of frag- 
mentation in the wave functions, hence collectivity. The 
precise knowledge of this residual interaction is therefore 
mandatory in mp — mh configuration mixing approaches, 
in particular when calculations are performed at spheric- 
ity. Consistently with the results for the wave functions, 
occupation probabilities display strong changes as can be 
seen in Table 1X1 in the case of proton and neutron 1(^5/2, 
2si/2 and 1(^3/2 orbitals. 



State Case 7rc?5/2 1^(1^/2 7rgi/2 7rc?3/2 z^c?3/2 

0+ (a) 5.053 5.047 0.603 0.606 0.344 0.347 

(b) 5.939 5.932 0.018 0.020 0.043 0.047 

0+ (a) 4.960 4.978 0.716 0.698 0.324 0.324 

(b^) 4.941 5.912 1.007 0.037 0.050 0.050 

(b^) 5.935 4.935 0.033 1.008 0.046 0.057 



TABLE X: Proton and neutron occupation probabilities of 
the (^5/2, Si/2 and ^3/2 orbitals, for 0^ and Oj states in ^^Si. 
Case (a) corresponds to a full calculation and case (b) to a 
calculation without residual proton-neutron interaction. 



IV. STATISTICAL PROPERTIES OF HIGHLY 
EXCITED CONFIGURATIONS 

One of the main issues raised by multiconfiguration ap- 
proaches is the number of relevant configurations for de- 
scribing low-lying states. Because of the proton-neutron 
excitations, this number can rapidly explode. In realistic 
calculations, one has to think about truncations based on 
physical arguments. When going beyond the mean- field 
with the mp — mh approach, it is assumed that short- 
range correlations have already been taken into account 
through the effectiveness of the nucleon-nucleon interac- 
tion used (the DIS Gogny interaction in our case). Our 
aim is to treat explicitly the long-range correlations cor- 
responding to the attractive part of the nucleon-nucleon 
interaction. Two standard types of truncation can be 
proposed, independently of the choice (or not) of a va- 
lence space: a truncation on the order of the excitations 
and/or a truncation on the configuration excitation ener- 
gies. Both types of truncations seem to be reasonable in 
the present context. In order to define appropriate trun- 
cations based on relevant physics argument, we discuss 
below the behavior of highly excited configurations. 

We have followed the direction used in Refs. [T3I - 
17] where the analysis was based on the properties of 
the strength function associated with the Slater deter- 
minants. Using second quantization and the standard 
Wick's theorem decomposition, the Hamiltonian Eq. 
([6]), is the sum of an independent particle part 1-Lq (one- 
body) and a residual part 1-L' (one-body and two-body), 

n = no^n'. (i4) 

The eigenfunctions \k) of the unperturbed Hamiltonian 

Mk)=ek\k), (15) 

describe noninteracting fermions. In the basis the 
residual interaction 1-L' has both diagonal, and off- 
diagonal, matrix elements: 1-L' = 1-L ^ 1-L . Full diag- 
onalization leads to the stationary states \a) and their 
energies So,-, 

n\a) =Sa\a). (16) 

The eigenfunctions \a) may have a complicated structure 
in the original basis 

\a)=J2^k\k). (17) 

k 

where A is a unitary matrix: 

^ A", A",' = 5"'^' , Yl ^kA",, = 5kk' ■ (18) 

k a 

A completely delocalized wave function \a) would have a 
number of relevant components close to the total dimen- 
sion N of the multiconfiguration space (for given exact 
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quantum numbers). In this limit the typical magnitude 
of each component is l/\fN. In general, a number No, 
of principal components \k) characterizes the delocaliza- 
tion of a state \a) in the given basis ([15]). Indeed, a 
two-body interaction can not couple configurations dif- 
fering by more than two particle states which implies a 
band-like Hamiltonian matrix, favoring the localization 
of eigenfunctions in the Hilbert space. Conversely, the 
fragmentation of simple basis states over the energy spec- 
trum can be provided by the strength function defined by 
the quantity 

Fu{E) = {k\5{E - H{p))\k) = \AtfS{E - £^). (19) 

a 

The strength function contains rich information but re- 
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FIG. 5: (Color online) Centroid energies (in MeV) for ^^Si. 
The X-axis corresponds to the number of centroids ordered by 
increasing energies. 

quires the knowledge of the full nuclear spectrum. Fortu- 
nately, one can study the main characteristics of the sys- 
tem from the first and the second moment of the strength 
function which does not require the actual diagonaliza- 
tion of H. These two moments are the centroid, E^, and 
dispersion, <j/c, of the state distribution, given by 

Ek= J dE EFk{E) = Ukk = efc + Uuk. (20) 

and 

al= fdE{E- EufFu{E) = ^ (n[,)\ (21) 

respectively. The centroid Ek coincides with the unper- 
turbed energy l-Lkk^ whereas the dispersion depends only 
on the off-diagonal matrix elements 



Below we discuss the example of the ^^Si isotope. Sim- 
ilar calculations have been done for other nuclei with sim- 
ilar conclusions. At this point, it is important to recall 
that the mp — mh configuration mixing formalism pre- 
sented in Section [IT] has been developed, in practice, using 
an axially-deformed harmonic oscillator basis. In order to 
preserve the spherical symmetry, all the mp — mh calcu- 
lations displayed in this study have been done for /3 = 0. 
In the following, our analysis is done in terms of different 
projections of a given angular momentum J. All con- 
figurations of the sd-shell (up to 12p — 12h) have been 
introduced in the wave functions; the maximum size of 
the Hamiltonian matrices that have been fully diagonal- 
ized was -90 000 x 90 000. 

Fig. [5] displays the values of the centroid energies Ej.^ 
where centroids are labeled by N/^ and ordered by increas- 
ing energy. One observes a characteristic behavior for all 
values of K from to 4. The lowest centroid energy, the 
one associated with the Op — O/i configuration, is -239.203 
MeV. Only few configurations have small excitation en- 
ergy. The level density increases rapidly with excitation 
energy, and most configurations are located in the range 
[-210; -160] MeV. The final increase, beyond -160 MeV, 
is an artifact of the finite valence space that reduces the 
number of possible highly excited states. 

Fig. [6] shows the dispersions Gk of Eq. ([2T]) for all con- 
figurations characterized by K = 0, 1, 2 and 3 and their 
excitation energies E* (in MeV). The color code indi- 
cates the number of configurations in a given excitation 
energy and (Jk bin. As expected, the structure of the 
configuration distribution is different for different values 
of K but most configurations are concentrated along a 
central line, and characteristic structures appear in the 
most dense areas. Such areas are more pronounced in 
the most dense zone in the K = case for all values of 
E*. A large majority of configurations display a disper- 
sion (Jk that is rather constant along the spectrum; its 
mean value can be evaluated to be — 9 MeV. This im- 
portant result which implies that strength functions can 
be characterized by essentially the same width is the first 
indication of expected exponential convergence. 

The strength function of ^^Si is presented in Fig. [7] for 
different values of iiT as a function of the energy counted 
from the corresponding centroid, E^ — E^. The value of 
the bin is set to 0.1 MeV. An average over all configu- 
rations has been done in order to reveal possible generic 
behavior. The upper panel displays the strength func- 
tion in linear scale and the lower panel in logarithmic 
scale. The central part of the distribution is intermedi- 
ate between a Gaussian and a Lorentzian one, whereas 
the behavior in the wings is found to be a decreasing ex- 
ponential. This is clearly visible in the logarithmic scale. 
In agreement with [13-17], the coupling of highly excited 
configurations with low-energy eigenfunctions therefore 
exhibits an exponential regime and the "3 a " rule char- 
acterizing the start of the exponential convergence, which 
is important in practical calculations of the level density 
by methods of statistical spectroscopy [3l|, seems to hold 
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FIG. 6: (Color online) Dispersions Gk (in MeV) of = 0+, 1^ 
the excitation energies E* of Slater determinants. The color 
excitation energy and dispersion. 

here also. 

In the case of the K = {) projection, one observes 
two unusual bumps at energies Eq^ — ~ 40 MeV 
and 60 MeV. Actually, the strength function (p!9]) com- 
bines two types of information: the density of config- 
urations and the mixing coefficients whose values are 
determined by the nucleon-nucleon interaction. In or- 
der to disentangle these effects, we calculate the den- 
sity of configurations p{E) = ^S{E — Ea) and the lo- 
cally normalized strength function {Fk{E)) / p{E) where 
{Fk{E)) = Y^f^Fk{E)/N. These quantities are shown in 
Fig. [8] (the density of configurations, upper panel (a), and 
the locally normalized strength function in logarithmic 
scale, lower panel (b)). The configuration density has a 
regular shape expected for a finite Hilbert space whereas 
the bumps of the strength function in Fig. [7] appear mag- 
nified on the normalized strength function plot (lower 
panel). One also notices that the coupling of matrix ele- 
ments is fully exponential (linear on logarithmic scale), if 
one forgets about the bumps. As the density of configu- 
ration is quasi-linear in logarithmic scale at large energy, 
one understands the exponential convergence behavior of 
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, 2^ and 3^ Slater determinants for ^^Si. The x-axis represents 
code indicates the number of Slater determinants per bins of 



the strength function observed on Fig. [71 The two bumps 
come from the specific interaction rather than from the 
statistics. The matrix elements that are responsible for 
the bumps correspond to the same proton-neutron ME 
discussed in section IIIIBI for ^^Si, which appear to be 
too large at the level of the approximate resolution of 
the mp — mh configuration mixing method presented in 
this study. Reducing by hand the values of these proton- 
neutron ME makes the bumps disappear. Hence, the 
study of the chaotic behaviour of highly excited Slater 
determinants can be considered as an additional way of 
highlighting not well-calibrated coupling ME. In other 
words, statistical properties may serve as an additional 
criterion in the validation process of phenomenological 
effective interactions. This result is consistent with and 
confirms our previous discussion of Section IIII B[ This 
analysis is similar in spirit to the statistical search of in- 
teraction matrix elements responsible for the equilibrium 
prolate deformation (32j . 

The exponential convergence observed in Fig. [8] is an 
interesting feature that might be evaluated analytically 
in the particular case of the finite-range Gogny interac- 
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FIG. 7: (Color online) The ^^Si strength functions associated 
with ii^^ = 0+, 1+, 2+, 3+ and 4+ components in linear scale 
(upper panel) and logarithmic scale (lower panel). Logarith- 
mic representation is shown in lower panel. Calculations have 
been performed with a 0.1 MeV bin. 



FIG. 8: (Color online) The ^^Si Slater determinant distri- 
butions p(E) (upper panel) and the ^^Si locally normalized 
strength functions {Fk{E)) / p{E) (lower panel) for = 
0^, 1^, 2^, 3^ and 4^ components . Calculations have been 
performed with a 0.1 MeV bin. 



tion. For example, taking the single-particle states as 
plane waves, the two-body matrix elements of the Brink- 
Boecker part of the Gogny force read: 

_ f d rid r2 g-ifci.rig-ifc2-r^g-(ri-r2)^//i^gifc3-rigifc4-r2 

J (2^)6 

_ r dVi ^-i(fei-fe2)-(ri-r2)^-(ri-r2)VM^ 
~ J (27r)6 ^ ^ 

. (22) 

where we have used the conservation law ki-\-k2 = ks-\-k4. 

Eq. (j22j) reduces to the Fourier transform of a Gaus- 
sian, 

(0102 |e-(---)^/^^ 103^4) = 4e-=^^, (23) 



with E = k'^ /2m and k = ki — k2 . Thus, the matrix ele- 
ment behaves, in the case of the Gogny interaction, as a 
decreasing exponential with respect to excitation energy. 
The value of the two ranges introduced in the Gogny in- 
teraction may serve as a guide to decide an upper limit 
where the exponential convergence regime settles. How- 
ever, this behavior exists for any physically reasonable 
interaction as revealed in shell model calculations 

Below we give three examples of exponential conver- 
gence behavior with increasing excitation energy of con- 
figurations, namely for the total energy, for the compo- 
nents of the correlated wave functions and for the occu- 
pation probabilities. Fig. [9] displays the evolution of the 
total energies of the ground and excited states in ^^Si, ac- 
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FIG. 9: (Color online) Convergence of ground state and a few excited state energies in ^^Si as a function of the number of 
configurations ordered by their centroid energies. 



cording to the number of configurations in the correlated 
wave function, ordered by increasing centroid energies. 
The four plots correspond to the different values of the 
total momentum projection K. In all cases, the total en- 
ergy changes rapidly when only a few Slater determinants 
with lowest energy are included. Then, for a larger num- 
ber of configurations a smooth regime settles, a behavior 
independent of the value of K. 

In Fig. dni the evolution of the global components, VF^, 
Eq. (fT3|) . of the wave functions of the 0]^ and 2^ states 
in ^^Si is displayed as a function of the number TVcent of 
centroids ordered by centroid energy. For n >9, the Wn 
are smaller than 10~^ and they are not shown. The evo- 
lution of the occupation probabilities for (^5/2, Si/2 and 
(i3/2 orbitals according to Accent is displayed on Fig. [TTl 
As can be seen from the two figures, both the Wn and 
the occupation probabilities display the exponential con- 
vergence. With all exponents being close, one can loosely 
interpret this behavior as a SI gnature of thermalization 
in a self-bound mesoscopic system. 



V. CONCLUSION AND PERSPECTIVES 

In this work, we have investigated the application of 
multiconfiguration methods for the description of low- 
energy nuclear spectroscopy. A few even-even silicon iso- 
topes have been studied, using the DIS Gogny interac- 
tion. At this stage, only the configuration mixing part of 
the method has been put into place; the renormalization 
of single-particle orbitals due to correlations has not been 
discussed in this paper. The results for the positive par- 
ity states in low- lying spectroscopy of ^^~^^Si have been 
found in rather good agreement with experiment, tak- 
ing into account the fact that the DIS Gogny interaction 
has not been a priori fitted to be employed in such a kind 
of approach. In particular, from the study of ^^Si, it has 
been found that, at the level of approximate resolution of 
the mp — mh configuration mixing method (no renormal- 
isation of orbitals) a few residual proton-neutron matrix 
elements of pairing type, that are not constrained in the 
fitting of the Gogny interaction, might disturb the repro- 
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FIG. 10: (Color online) The W^ components (Eq. {[T3))) of 0^ and 2+ wave functions in ^^Si. 



duction of excitation energies. In relation to this, the 
importance and magnitude of the proton-neutron resid- 
ual interaction has been discussed. 

The question of the pertinent configurations that have 
to be introduced in the mixing has also been addressed. 
Statistical generic behavior of highly excited configura- 
tions has been put forward. In particular, the exponential 
convergence already revealed in shell model studies has 
been confirmed in our approach. This is an encouraging 
feature that may help to handle the very large number of 
configurations that appear in multiconfiguration meth- 
ods, in particular in nuclear physics where two kinds of 
particles exist. 

The present study proposes interesting and challenging 
issues at different levels. The renormalization of orbit als 
under the influence of correlations is an important ques- 
tion that will be analyzed in further studies. In atomic 
physics, this renormalization has been proved to play a 
key role in strongly correlated systems. Even though the 
associated orbital equation, see Eq. ([8]), looks simple, its 
exact solution is far from being trivial. At present, in the 



most advanced applications to atomic physics, it is solved 
approximately as the correlation term that depends on 
the two-body density matrix can be very complicated. 

A second issue concerns the improvement of the Gogny 
interaction in order to be able to use it not only in HFB, 
RPA and GCM-type methods but also in multiconfigu- 
ration approaches. Work is in progress in this direction. 

The third issue deals with the generic behavior of 
highly excited configurations. The exponential conver- 
gence and corresponding extrapolations can be of con- 
siderable help for controlling in a safe way possible new 
truncation schemes introducing explicitly only pertinent 
configurations. To this aim, a formalism of the Feshbach 
type projection might be quite useful. It is worth men- 
tioning that the analysis of statistical regularities allowed 
us to identify specific matrix elements of the Hamiltonian 
responsible for spectroscopic inadequacy (discussion on 
30Si). 

For the specific goal of nuclear spectroscopy, the inves- 
tigation of transition probabilities and negative parity 
states in silicon isotopes would be of great interest and 
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are fully occupied for each isospin; 

- a set of active orbitals containing — particles for 
each isospin; 

- a set of unoccupied higher-energy individual states for 
each isospin. 

Under this prescription, the proton and neutron Slater 
determinants, Eq. ([3j), are defined as 



\^-r)= n 4ll4\0)= n (Al) 

where \^c-r) is the part of the wave function that describes 
the core. 

In addition, the Hamiltonian H{p)^ Eq. ([1]), can be 
written by separating isospin contributions as 



H{p) = J2 H- + V^-{p). 



(A2) 



FIG. 11: (Color online) Evolution of the occupation proba- 



bilities of proton 



and 



^'3/2 



orbitals as function of 



Then, from Eq. (|6|), 

= E + E ^^mnatmarn] + (A3) 



may provide essential information on the residual inter- 
action and properties of single-particle orbitals, including 
nuclei far from stability. 
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Appendix A: "Core + valence space" formulation 

Eqs. (j2j) and (j5j) are quite general. As mentioned in the 
Introduction, the Hilbert space has to be truncated in any 
realistic calculation. Several truncation schemes, sup- 
ported by physics arguments, can be utilized, for example 
a limitation on the excitation order of mp—mh configura- 
tions, a limitation in the number of single-particle states 
used for the configuration mixing, etc. This appendix 
deals with a truncation that corresponds to the descrip- 
tion in terms of "core + valence space" . In this approach, 
the system, that comprises nucleons of each isospin, 
is separated into: 

- an even-even core where the lower individual states 



where the generalized rearrangement coefficients TZJj are 
given by 

T^Ln = J rrm{r,cT)<Pru{r,a)B{r)d^r, (A4) 

with (j)ri the single-particle wave functions. The difficul- 
ties in the calculation of 7^^^ come from the evaluation 
of the rearrangement field B{r) from the correlated wave 
function ([2]): 



B{r) = E^<<^«^«^ E^^<^< 



ijkl 



dVipl 
' dp{r) 



I ^CKtt ^Oiu ) 



(A5) 

The expressions (|A2[) . (|A3[) and (|A5[) are specific to the 
DIS Gogny interaction as only the proton- neutron terms 
are generated by the density-dependent part of the inter- 
action and the density-dependence is a contact interac- 
tion. In the limit of no configuration mixing, one recovers 
the standard HE expression for the rearrangement term. 

Erom Eq. (|A3[) . one derives explicit expressions for 
^a-,,a^,a'^a'^ in tcrms of the "core + valence space" for- 
mulation. In the following evaluation oH-La^^au^cx' a' , only 
the terms carrying the core contributions are given. 



1. Proton and neutron diagonal contributions 

The proton and neutron diagonal contributions are ob- 
tained for l^a^) = l^a;) and \^a^) = \^a'^) as 



n 



(A6) 
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which yields 



1 



(A7) 



In Eq. (|A7[) . Zq,^ stands for the occupied single-particle 
state in the \^a-r) Slater determinant. The "core + va- 
lence space" separationthen leads to: 



(A8) 



In Eq. ()A8p . the first and third terms are pure core 
contributions. The second and forth terms are the con- 
tributions from the valence space. The mixed fifth term 
includes single-particle orbitals of both the core and the 
valence space. 

2. Proton-neutron diagonal contribution 



In Eq. ()Alip , the first term comes from the core and the 
second one from the valence space. The last two terms 
express the coupling between the core and the valence 
particles. 



3. Proton and neutron non-diagonal one-body 
contributions 



The proton or neutron one-body contributions arise 
when \^a^) and 1^^;) differ in one particle state and 
\^a^) = \^a' ) or when \^a^) and ) differ from one 
particle state and = 1^^' ), respectively. Then, the 

contributions to H are 

ATJ 



(A12) 



,=Ar-+i 



and 



AT- 



E {<j<\v{p)\ 



(A13) 



.=Ar-+i 



In Eqs. ()A12p - (IA13p . the indices j and / refer to the 
single-particle state different between the configurations 
a and a' . The first term in Eqs. ()A12p - (IA13p corre- 
sponds to the core contribution and the second one to 
the contribution of the valence space. 



The proton-neutron diagonal contribution is obtained 
for |<|)c,;) and |<I>«J |<1><), as 

K:«.,c..«. = {^aMVn^a.^a.). (A9) 
Expansion of Eq. ()A9p yields 

AT^ AT^ 

k:c..,c..«. = Y.Y.^'-jaAv\CX.)- (Alo) 

Making the "core + valence space" separation leads to 

AT- AT^r 

AT^ AT^ 

+ E E (ic.Jo.AV\ 

=Ar- + l =N^ + 1 /Am 
AT- i^-L-Lj 

AT- AT- 
+ E E (^«.i«JV'|z«.jc,J. 



4. Proton-neutron non-diagonal one-body 
contribution 

As for the proton and neutron one-body contribu- 
tions, the proton-neutron one-body contribution arises 
when l^a-j,) and ) differ by one particle state and 
|<I>Q,^) = \^ai) OT when l^a^) and \^ai) differ by one par- 
ticle state and \^a^) = Similarly to Eqs. ()A12p - 
()A13p . the proton- neutron non-diagonal one-body contri- 
bution is then given by 

^--^ (A14) 

+ E {^c.JaJV^''{p)\^^A' 

i«^=Arj+i 

In Eq. (|A14p . the indices j and / refer to the single- 
particle states that differ between the configurations a 
and a' . Still, the first term describes the core contribu- 
tion and the second one the contribution of the valence 
space. 
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